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We study the inverse boundary crossing problem for diffusions. 
Given a diffusion process Xt, and a survival distribution p on [0, oo), 
we demonstrate that there exists a boundary b{t) such that p{t) = 
P[r > t], where r is the first hitting time of Xt to the boundary b{t). 
The approach taken is analytic, based on solving a parabolic varia- 
tional inequality to find b. Existence and uniqueness of the solution to 
this variational inequality were proven in earlier work. In this paper, 
we demonstrate that the resulting boundary b does indeed have p as 
its boundary crossing distribution. Since little is known regarding the 
regularity of 6 arising from the variational inequality, this requires a 
detailed study of the problem of computing the boundary crossing 
distribution of Xt to a rough boundary. Results regarding the formu- 
lation of this problem in terms of weak solutions to the corresponding 
Kolmogorov forward equation are presented. 

1. Introduction. Let {Bt}t>o be a standard Brownian motion defined 
on a filtered probabifity space (r2,P, {J-^t}t>o) satisfying the usual conditions. 
We consider a diffusion process {Xt}t>o defined by the stochastic differential 
equation 

dXt = fi{Xt,t)dt + a{Xt,t)dBt Vt > 0, 

where /i, a : M x — )■ M are smooth bounded functions, with bounded deriva- 
tives^ and infigxK+ cr > 0. We assume that Xq is independent of B, and has 
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The main results in the paper are true under less restrictive assumptions on the 
coefficients, although we have not sought to determine the minimal assumptions under 
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initial distribution P(Xo < = po{x, 0), with density po{x, 0). For y € M and 
t > s, we further denote by F{y,s;-,t) and p{y,s;-,t) the transition distri- 
bution and density, respectively, of Xt given Xg = y: 

F{y, s; x,t):= F{Xt <x\Xs = y), p{y, s; x, t) := ^^^^^^'"^'^^ 

In the sequel, we denote by po{-,t) and po{-,t) the density^ and cumulative 
distribution of Xt : 

Po{x,t):=F{Xt<x)= [ po{y,0)F{y,0;x,t)dy, 

JR 

po{x,t):= — := / Po{y,0)p{y,0;x,t)dy. 

For a given function 6:M+ — t- [—00,00), the first boundary crossing time 
f is defined to be 

(1.1) f = mi{t>0\Xt<b{t)}. 

We shall also have occasion to consider the related, but less commonly used, 
first time that Xt goes strictly below b: 

(1.2) T = mi{t>0\Xt<b{t)}. 
We are interested in the following two problems. 

1. The boundary crossing problem: for a given function 6:R_|_ — )■ [—00,00), 
compute the survival distribution of the first time that X crosses b; that 
is, 

(1.3) p{t)=F{f>t). 

In this case, we denote p = V[b]. 

2. The inverse boundary crossing problem: for a given survival distribution 
p on (0,00) find a function b such that b satisfies (1.1), (1-3). If such a 
boundary exists and is unique, we denote it by 13\p]. 

The boundary crossing problem is classical, and the subject of a large liter- 
ature. The inverse boundary crossing problem has recently been the subject 
of increased interest by probabilists and researchers in mathematical finance. 
The main purpose of this paper is to show that the inverse boundary crossing 
problem is well posed. 



which they hold. The stated assumptions ensure that the apphcations of various results 
from the theory of PDEs, a change of variables, and an application of Ito's lemma are all 
valid. 

^At t — 0, po is interpreted as the distributional derivative of po- 
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According to Zucca and Sacerdote (2009), the problem was originally 
posed by A. N. Sliiryayev in 1976, for the special case where Xt is a Brow- 
nian motion and p is the exponential distribution. Dudley and Gutmann 
(1977) and Anulova (1980) showed that there exists a stopping time with 
the given distribution; however, this stopping time is not realized as the 
first time the process X crosses a boundary h. Recently, there has been an 
increase of interest in the problem due to its importance in applications. In 
mathematical finance, with Xt an indicator of a firm's financial health, and 
p the distribution of its time to default (estimated from the prices of market 
instruments), the problem is to find a default barrier that reproduces the 
given default distribution. Many authors have proposed numerical methods 
for finding such a boundary, including Hull and White (2001), Iscoe and 
Kreinin (2002), Huang and Tian (2006), Avehaneda and Zhu (2001) and 
Zucca and Sacerdote (2009). A formulation of the problem in terms of non- 
linear Volterra integral equations has been given by Peskir (2002) [see also 
Peskir and Shiryayev (2005) for a more detailed discussion]. 

The numerical method proposed by Avellaneda and Zhu (2001) is most 
relevant to our work. They note that for sufficiently smooth boundaries 6, 
the function U{x,t) = 9xP(t >t,Xt< x) should be the solution of the free 
boundary problem 

{CiU{x,t)=0, for x>6(t), t>0, 

U{x,t)=0, for x<bit), t>0, 

U{x,0)=po{x,0), forxGM, 

with the free boundary condition 

(1-5) P{t) = -'^{a^UU,=b{t) Vi>0, 

where Ci is the differential operator 

Avellaneda and Zhu (2001) perform a change of variables to "straighten 
out the boundary," and then solve the resulting transformed PDE numeri- 
cally using finite differences. 

An analytic study of the inverse boundary crossing problem was initiated 
by Cheng et al. (2006). In that work, we defined 

(1.7) w{x,t)= / Uiy,t)dy. 

J X 

Formally, direct calculation from (1.4), (1.5) shows that w should satisfy the 
free boundary problem 

(Cw{x,t) = 0, for X > 6(t), t > 0, 

w{x, t) = pit), for X < b{t), t>0, 

Wxix,t)=Q, for X < 6(t), t > 0, 

. w{x, 0) = 1 — po{x, 0), for X G M, 



(1- 
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where C, is the differential operator 
Formahy, 

(1.10) w{x,t)=¥{f>t,Xt>x), 

providing the connection between the probabihstic problems (1.1), (1-3) and 
our analytic approach. Based on the free boundary problem (1.8), one can 
infer that w should satisfy the variational inequality 

m&x{Cw,w-p} = in L°°(]R X (0,c5o)), 

(1.11) 

w{;0) = l-poi;0) onM, 
and b can be recovered from w by 

(1.12) b{t):=mi{x\p{t)-w{x,t)>0} Vt > 0. 

In Cheng et al. (2006), the existence and uniqueness of a viscosity solution 
to (1.11) was proved. However, no attempt was made to connect the re- 
sulting functions w,b to the original probabilistic formulation of the inverse 
boundary crossing problem. In this paper, we show that b does in fact give 
a boundary that reproduces the survival distribution p. This is complicated 
by the fact that it is very difficult to prove the regularity of the boundary b 
derived from the variational inequality.^ As a consequence, in order to verify 
that b has the required hitting distribution, we must first study the problem 
of computing the boundary crossing probabilities of diffusions to nonsmooth 
boundaries. To this end, for a given function b : (0, oo) — )• [— oo, oo), we define 

(1.13) b*(t):=max\b(t),]hRb(s)}, b*^(t) :=]hRb(s) Vt > 0. 

I s-i-t J s yt 

When needed, we also define b*{0) := lim^x^o ^(•s)- We also employ the nota- 
tion Qf, := {{x,t)\x > b{t),t > 0}. It will turn out that the inverse boundary 
crossing problem is most naturally formulated in the following spaces: 

So:=|fe:(0,oo)^[-oo,oo)|6=6*=6*„,p('|J f] {X ^ > b{s)}] = l\ , 

Po := {p £ C{[0, oo))|p(0) = 1 > p{s) > p{t) >Oyt>s> 0}. 
The main result of this paper is the following theorem. 



®The problem of the regularity of the boundary has subsequently been investigated by 
Chen (2011). 
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Theorem 1. 

1. For every p ^ Pq, there exists a unique viscosity solution, w, for the sur- 
vival distribution of the inverse boundary crossing problem associated with 
p [i.e., a viscosity solution of problem (1.11)]. In addition the unique so- 
lution, w, satisfies 

< 1 - [w{x,t) +po{x,t)] < l-p{t), 

(1.14) 

w{x, t) < w{y, t) Vt > 0, X G M, X > y. 
Consequently, the operator B 
(1.15) b{t) = B\p]{t) := inf {x e R\w {x, t) <p{t)} Vt > 0, 

is well defined on Pq. 

2. For every p G Pq, B[p] € Bq and {roB)[p] =p, where {V[h]){t) =¥{t > t) 
and f is defined as in (1.1). 

(This implies that for a given p G Pqj ^ '■= ^[p\ is a solution of the 
inverse problem since V\f)\ = {V o B)\j)\ = p.) 

3. For every b G Sq, V[b] G Pq and {B o 'P)[bj = b. 

(This implies that for a given p Pq, if b & Bq is a solution of the inverse 
problem, i.e., V[b] = p, then b can be identified as b = {B o V)[b] = B[p], 
the viscosity solution of the inverse problem.) 

4. // {w,b) is a classical {i.e., w + pq G C(M x [0,oo)), dxvu G C(M x (0,cx3), 
dtw,dxxW G C{Qf))) solution of the free boundary problem (1.8), then b is 
the solution of the inverse boundary crossing problem associated with p, 
that is, ^[6] =p. Similarly, if {U,b) is a classical {i.e., U — pQ € C(M x 
[0,cxd)), dtU,dxxU GC{Qb)) solution of the free boundary problem (1.4), 
(1.5), then V[b] =p. 

The proof of the above theorem proceeds as follows. We begin by study- 
ing the direct problem of computing the distribution of f, and the func- 
tion w{x,t) = P(f >t,Xt> x) for boundaries b G Bq. By considering a care- 
fully constructed discrete approximation scheme motivated by (1.4) when 
b is known, we are able to show that w is the unique viscosity solution to 
(1.8). Elementary calculations verify that the viscosity solution of the varia- 
tional inequality (1.11) also solves (1.8). Once we have also determined that 
{x > b} = {w < p}, the verification proceeds by relatively straightforward 
arguments. 

We note that the sequence of stopping times constructed in our discrete 
time approximation actually converges to the first time that X is strictly 
below the boundary b, t as given by (1.2). This definition of the boundary 
crossing time is slightly different from the standard one (1.1) for f. We 
have found that for the analytic approach we take here (particularly for 
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rough boundaries), our definition is more convenient. In Section 2 below, 
we show that for boundaries with minimal regularity properties (including 
those arising from the solution to the inverse boundary crossing problem, 
be Bo), P(r = f) = l. 

The remainder of the paper is structured as follows. The second section 
proves measurability properties of r and f, and proves that these times are 
almost surely equal. In addition, it gives preliminary results that are needed 
for the study of our approximation scheme. The third section studies the 
approximation scheme in detail, and proves convergence. The convergence 
provides a rigorous connection between the probabilistic definition of the 
survival probability p and the PDE definition of the survival distribution w. 
The fourth section formulates viscosity solutions for the direct problem of 
computing p for a given 6 G -Bq , and shows that the survival distribution w 
gives the unique viscosity solution for the direct problem. The fifth section 
provides the link between the variational inequality studied in Cheng et 
al. (2006) and the inverse boundary crossing problem. It also provides a 
sufficient condition under which the resulting boundary b is continuous. 

2. Crossing times of upper-semi-continuous boundaries. We calculate 
boundary crossing distributions for rough boundaries based on discrete time 
approximations to be studied in the next section. In order to ensure conver- 
gence of our approximation scheme, the time points used must be chosen 
carefully. We refer to the selected points as the "landmark points" of the 
boundary. In this section, we begin by defining the landmark points and 
investigating their properties. Using these properties, we study the measur- 
ability of r and f, show that the boundary crossing times of b and b* are 
equal and that P(r = f) = 1 for b G Bq. 

Definition 1. Let 6: (0, oo) [—cx), oo), and 6* be its upper-semi-continuous 
envelope. The set of landmark points of b, denoted by A(6), is defined as 
follows: 

(2.1) A(6) := U An{b), A„(6) := {ti\i G N}, 

nSN 

(2.2) i^=inf|te 

The following lemma summarizes some properties of the landmark points 
that are used in the paper. 

Lemma 2.1. Let b: (0, oo) — )■ [—00,00), and let its landmark points A(6) 
be defined as in (2.2). 



i + l 



b*{t) > 



sup 



b{- 
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1. For i,nGN, b*{ti)>b*{s) for every s € [2~"i, 2-"(i + 1)). 

2. For i,n£N, either = t^^^ or = so A„(6) C A„+i(6). 

Proof. 1. If s G (2""i,2-"(i + 1)), then by definition 6* (4) > b{t) for 
t € (s — e, s + e) with e small enough. > b*{s) follows immediately. If 

s = 2-"i, and b*{s) > > sup{6(t)|t G [2-"z, 2-'^(i + 1)]} then s = ^ by 

(2.2). 

2. We first claim that if G [2-"2, 2-("+i)(2z + 1)) then tl, = tl\^. Clearly, 
> sup{6(s)|s G [2-"i,2-"(i + 1)]} > sup{6(s)|s G [2-"i, 2-("+i)(2i + 
1)]}, so by definition < t^. If the inequality is strict, there is a (5 > small 
enough so that + C (2-"z,2-("+^)(2i + l)), and since b*{tl\^)>b 

on this interval, we obtain — ^*(*n)' contradicting the definition 

of 4. A similar proof shows that if tj, G [2-('^+^)(2i + l),2-'^(i + 1)), then 
tn = tn+i- Finally, it is easy to see that if tjj = 2~"(« + 1), then sup{6(s)|s G 
[2-"i,2-"(i + l)]} = sup{6(s)|s G [2-("+i)(2i + l),2-"(i + l)]} > sup{6(s)|s G 
[2~"'i, 2~(""*"^)(2z + 1)]}, after which repeating the same argument by contra- 
diction ensures that t^^ = t'^^^ ■ □ 

The following lemma collects some properties of upper-semi-continuous 
functions that are used throughout the paper. The proofs are elementary, 
and are omitted. 

Lemma 2.2. Let b: (0,oo) — )• [—00,00) be upper- semi- continuous. 

1. If x: [0, 00) (—CO, 00) is continuous, then for all t>0, 

inf{s > 0|a;(s) < b{s)} > t <^ x{s) > b{s) Vs G (0, t] 

2. The set Q;, := {{x,t)\x > b{t),t > 0} is open. 

The next proposition addresses two main issues. First, it considers the 
measur ability of r and f, to ensure that the various functions considered in 
the remainder of the paper are well defined. Second, it shows that for the 
purposes of computing the distribution of r, it is enough to consider the 
upper-semi-continuous envelope, b* , of the boundary b. We observe that the 
result for r holds with minimal assumptions on the function b (we have not 
even assumed measur ability) . Furthermore, we note that some of the results 
on measurability could be derived by applying more general theorems [e.g., 
f is the first hitting time of the two-dimensional process {Xt,t) to the set 
{{x, s)\x < b{s)}, which is closed when b is upper-semi-continuous]. However, 
we have decided to present elementary proofs of these assertions to make 
the paper more self-contained. 
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Proposition!. Let b:{0,oo) ^ [-00,00). 

1. Let b* be as in (1-13) and A{b) be as in (2.1) and (2.2). Then for every 
t>0, 

(2.3) n {X,>bis)}= n {Xs>b*is)}. 

sG{0,i) se(o,t)nA{6) 

Consequently, we can define the first boundary crossing time r : O — )• 
[0,00], the survival probability p: [0, 00) — > [0,1], and the survival distri- 
bution w.Rx [0, 00) — > [0, 1] by 

(2.4) T{uj):=mf{s>0\Xs{uj)<b{s)} ^uGQ, 

(2.5) p{t) ■.= F{Xs>b{s) yse{0,t)) Vt>0, 

(2.6) w{x,t) ■.= F{Xs>b{s) yseiO,t),Xt> x) VxG]R,t>0. 

In addition, r is an optional time with respect to the filtration generated 
by the process X, {r > t} € Ti^yt > 0. Also 

p{t) = F{T>t), w{x,t)=F{T>t,Xt>x). 

2. Let the (conventional) first crossing time f:Q ^ [0, 00] be defined by 

f{uj):=mf{s>0\Xsiuj)<b{s)} Vo; € Jl. 

If b is upper- semi- continuous, i.e., b = b*, then t is a stopping time with 
respect to the filtration generated by the process X , {u ^ Q\f{uj) > t} (z 
yt > 0, so that we can define p{t) := P(f(w) >t)yt> 0. 



Proof. 1. When t = 0,we have {r > 0} = J] = G ^\Xs{uj) > 6(s), Con 
Vs G (0,t)} and p{0) = F{Q) = 1. Now we assume that t > 0. It is easy to 
verify that 

n > Ks)} = {Xs > bis) Vs € (0, t)} = {T> t}. 

sG(0,t) 

Hence, to complete the proof of the first assertion, it suffices to verify (2.3). 
Suppose uje{Xs> b{s) Vs G {0,t)}. Then Xs{u)) > 6(s) for every s G (0,t). 
For every .s G (0,t), by the continuity of X,{uj), 

Xg{uj) = lim Xs{uj) > max "I b{s), lim b{s) > = b*{s). 

As s G (0,t) is arbitrary, we have lo G {Xg >b* Vs G (0,t)}. Thus, {Xg > 
6(s) Vs G (0,t)} C {Xs > b*{s) Vs G (0,t)} C {X, > b*{s) Vs G (0,t) n A(6)}. 
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Next, suppose oj G > h*{s) Vs G (0,t) n A(6)}. Let s G (0,t) be 

arbitrary. We want to show that Xg > b{s). For each integer n satisfying 
2"" < s, let in be the integer such that s G [i„,2~"', (i„ + 1)2""). Then 4" G 
A(6) and ^^^^1^) > 6*(4") > b*{s) by Lemma 2.1. Hence, 

Xsiuj) = hm Xi„(w) > lE^ 6* (4") > r(s) > b{s). 

n— >-oo n— >oo 

Since s is arbitrary, we see that w G {X^ > b{s) Vs G (0,t)}. Consequently, 

{r>t}= n {x,>6(s)}= n {^.>^*(^)} 

se(o,i) se(o,t) 
= n {Xs>b*{s)}€Ti'. 

s€{0,t)nA(b) 

This proves (2.3) and also the first assertion. 

2. Assume that b is use, that is, b = b* . Then by the continuity of the 
sample paths of X, if s > and Xg{uj) > 6(s) then there exists 6 > such 
that Xs{uj) > b{s) + 6 for all s G [s — (5, s + 5]. By the Heine-Borel theorem, 
if Xg > b{s) for every s G [a,c] C (0,oo), then there exists a large integer i 
such that Xs > b{s) + 2~* for every s G [a, c]. Hence, for every t > 0, 

{f>t}= n {^s>K«)}=n n {^s>b{s)} 
(2-7) =nu n {Xs>b{s)+2-^} 

ne'MieN se[2-"t,t] 

= n {X,>6(.) + 2-}GJ-f. 

ngNigNsG{2-"t,t}U(A(6)n[2-"t,t]) 

This completes the proof of the second assertion. □ 

The following proposition justifies our choice to work with r, the first 
time the process is strictly below the boundary, rather than f, the first time 
the process hits the boundary. In particular, the second assertion implies 
that these times are almost surely equal, and hence they have the same 
distributions (so solving the inverse boundary crossing problem for r is the 
same as solving it for f). We will see in the next section that it is easier to 
work with r in defining approximations for the boundary crossing problem. 

Proposition 2. Letb: (0, oo) [-oo, oo) . 
1. Let b* and bt be as in (1.13). Then 

{w G Q\t{uj) = t or f{uj) =t}c{uoe n\Xt{Lo) G [b*_{t),b*{t)]} Vt > 0. 
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2. The set {r 7^ f} has probability zero, so that 

p{t)=p{t-)=p{t-)=F{T>t) yt>o, 

p(t)=p{t+)=p{t+)=F{T>t) \ft>0. 

Consequently, if b* =b*_, then p G C{{0, 00)), p = p on (0,oo), and p € 
C([0,oo)). 

Proof. To prove the first assertion, let t > and w G {r = t} U {f = t}. 
Then Xs{io) > b{s) for all s G (0,i) so Xt{uj) > Ihil^/.t ^(s) = b*_{t). Also, 
t{uj) = t 01 f{u)) = t implies that there exists a sequence {si} of positive 
numbers such that linii^oo Si = t and Xs^ {oo) < b{si) for all i. This im- 
plies that Xt{uj) = liuii^oo Xs^{uj) < limi^oob{si) < b*{t). Hence, Xt{uj) G 
[b*_{t),b*{t)]. Also, note that if b*_{t) = b*{t), then P({r = t} U {f = t}) < 
P{{Xt = b*it)}) = 0. 

Since the family {t > t + e}e>~t of sets is monotonic in e, we see that 

p{t+) = lim p{t + e) = lim P(t >t + e) 

e\0 e\0 



|J{r >t + e} j =P(r>t) Vt > 0, 

e>0 ^ 



p{t-) = lim p{t -e) = lim P(t > t - e) 

£\0 e\0 



(f]{T>t-£}\=F{T>t)=p{t) Vt>0. 



V>0 

Clearly, if b*^{t) = b*{t), then p{t-) -p{t+) = P(r = i) = so p is contin- 
uous at t. Similarly, when b* = b so p \s well defined, we have p{t) = P(f > 
t) = p{t+) and p{t-) = P(f > t). 

To complete the proof, it remains to show that the set {r 7^ f } has prob- 
ability zero. For every w G 

t{uj) = inf{s > 01X3(0;) < b{s)} > inf{s > 01X^(0;) < b{s)} = f{uj). 

Now, suppose t{ijj) 7^ ^^(a;). Then we must have f(u;) < t{uj). Set t = t{uj). 
Then Xg > b{s) Vs G (0, t). By continuity, we also have Xs > b*{s) Vs G (0, t). 
Set t = t{u}). If t = 0, then by definition, there exists r G (0,t) such that 
Xr{uj) < b{r). If t > 0, then by definition, there exists r G [t,t) such that 
Xr{io) < b{r). As Xs{uj) > b*{s) for all s G (0,i), in either case, we have 
r G (0,t) and Xriuj) = b*{r). Taking ri G A(6) n (0,r) and r2 G A(6) n (r,t) 
we obtain minsgjri,r2]{^s('^) — ^* (•?)} = 0. Hence, 

{r/f} = {f<r}c U U B(ri,r2), 

rieA(b) r2eA(6)n(r-i,oo) 
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where for every < a < c < oo, 

B(a,c) = lujen\ mill {^^(a;) - b*{s)} = o). 

s&[a,c] ) 

Note that for each c> a > 0, B(a,c) is J-c measurable since 
B(a, c) = ( mill {X, - fe*} > o) 

oo 

\ U{min{X,-6*(s)}>2-"|, 

\ , Ug a.c J 

n=l 

f min{X, -6*(s)}>/i|= H {X, > 6*(s) + /i} G V/i G M. 

•-sGfa.cl J ' ' 

s£{a,c}U{A{b)n[a,c]) 

It remains to show that for each c > a > 0, the set B(a,c) has measure 
zero. Suppose, on the contrary, that P(B(a,c)) > for some fixed c> a > 0. 
Fixing to £ (0,a), we then have 

0<P(B(a,c))= / P(B(a,c)|Xto =z)po(^,to)d^- 
Consequently, there exists a finite number M > such that 
/ F{B{a,c)\Xt,=z)po{z,to)dz>0. 

For each /i € M, we consider the set 

B'*(o,c) = | min (X, -6*(s)) = /ij. 

l-sGla.c] J 

For the process {Xt}t>to, for each uj G {X^ = z} and /i G M, we denote by 
uj'^ the element in {Xtg = z + h} such that Xt{uj^) = h + Xt{ijj) Vt G [to, oo). 
Then 

min -6*(s)) = ^ min (X,(a;'') - = /i. 

s6[a,c] sG[a,c] 

Assume for simplicity that we are dealing with Brownian motion. [By a 
change of variables, we can assume £7 = 1; see Section 4. Since fi is smooth 
and bounded, if {Xt} is not a Brownian motion, we can use the Girsanov 
theorem [Karatzas and Shreve (1996)] to change to an equivalent measure 
under which Xt is a Brownian motion, and the argument below can still be 
used to show that P(B(a, c)) = 0.] By the translation invariance of Brownian 
motion and the Markov property, we have 

P(B(a,c)|Xt„ =z)=P(B'^(a,c)|Xt, =z + /i). 
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Hence, 

F{B^{a,c))= [ ¥{B''{a,c)\Xt,=z + h)po{z + h,to)dz 
= [ ¥{B{a,c)\Xt^=z)poiz + h,to)dz 

JR 

> min P^^^ph^ r F{Bia,c)\Xt,=z)poiz,to)dz>0. 

Note that all elements in {B^{a,c)}fi£R are disjoint and measurable. We 
then obtain a contradiction since does not contain an uncountable disjoint 
union of measurable sets with positive probability. Thus, B(a,c) must have 
probability zero, for every pair (o,c) with a > c > 0. Consequently, the set 
{f 7^ r} has probability zero. This completes the proof of Proposition 2. □ 

Theorem 2. The operator V defined by 'P[b]{t) = P(r > t) maps Bq 
toPo. 

Proof. Suppose b € Bq. Then b* = b*_ sop := V[b] G C((0,(X))). In ad- 
dition. 



lim^)(e) = limpf f] {X, > 6(s)} ) 

U n {Xs>b{s)})=l=p{0). 



£>Ose(0,e) 

Hence, p € C([0, oo)). It remains to show that p > on [0, oo). Since 
p(0+) = p(0) = 1, there exists e > such that p{t) > for every t € [0,e]. 
Let T > e. The upper-semi-continuity of b implies that M := sup^gjo.T] ^(^) 
is finite. Then P(r > e,X^ > M) > 0. Using standard results for a constant 
barrier 6 = M on the set {r > e, > M} for the time interval [e, T] , we see 
that P({r >e,Xs> M yse [e, T]}) > 0. Hence, p{T) > 0. As T is arbitrary, 
we see that p > on [0, oo), so that p £ Pq. □ 

Proposition 3 (A semi-continuous dependence property of V). Assume 
that 6, 6i, 62, . . . are upper- semi- continuous functions having the property 

fei<&2<&3<---, b{t) = Urn bn{t) yt>0. 

n— >oo 

Let p = V[b] and pn = P[6„]. Then for every t > 0, p{t) = lim„_j.ooPn(i)- 
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Proof. For every t>0, 

{r>t}= fl {Xs>b{s)}= fl f]{Xs>bn(.s)} 

se(0,t) s6(0,t)nGN 

= n n {Xs>bn{s)}. 

neNsG(0,f) 

Hence, 




= lim Pn{t). □ 

n— >oo 

3. Approximating sequences for boundary crossing times. In this sec- 
tion, we use the landmark points to construct straightforward approxima- 
tions that eventually, upon passing to the limit, will allow us to transfer the 
problem of calculating the survival probability to problems of solving partial 
differential equations. The advantage of studying the first time the process is 
strictly below the boundary is suggested by comparing the relative complex- 
ity of the expressions (2.3) and (2.7). In this case, a simple approximation to 
the survival probability and distribution can be developed. We approximate 
a real barrier 6 by a simple barrier bn defined by bn{t) = b*{t) if t G A„(6) 
and bn{t) = — oo otherwise. The approximate problem then involves only the 
random variables {Xt\t G A„(6)} so that all relevant probabilities can be cal- 
culated through transition probability densities. Though it turns out that 
survival probabilities computed using both viewpoints are equivalent, we do 
not see a simple adaptation of our method that allows us to approximate f 
directly without appealing to the results in Section 2. 

Proposition 4. Let b: [0,oo) — ?> [—00,00) be upper- semi- continuous and 
{t,p,w) be defined as in (2.4)-(2.6). Let A{b) = Un6N''^"(^) landmark 
points of b, An{b) = G N} and 

(3.1) Tn{u;):=mm{seAn{b)\Xs{u;)<b{s)}, Pn{t) ■.= r{Tn>t), 

(3.2) Wn{x,t):=F{Tn>t,Xt>x). 

Then the following hold: 
1. For all (x,t), 

/•oo 

Wn{x,t) = / Un{y,t)dy, 

J X 
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where, for t S [0,t°],{r„ > t} = Vl,Pn{t) = 1, -(/;„(•,*) = 1 -po(-,i) and 
Un{-,t) = po{-,t), and when te (t^,*^"^^] withk^n, {rn > t} = fljLil-'^t^ > 
b{tU} and 

POO 

Un{x,t)= / Uniy,t'^)p{y,t^;x,t)dy, 

(3.3) 

^ ' p poo 

Pn{t)= / Un{y,t)dy= / U{y,ti)dy. 

2. For every n G N, T„ > Tn+i >T,Pn> Pn+l >P,Wn> Wn+l >'W,po>Un> 
Un+1>0. 

3. There exists [/ : M x (0, oo) [0, oo) suc/i that for every co £ Q,,t > and 

lim {Tniu}),Pn{t),Wn{x,t),Un{x,t)) = {t{uj) , p{t) , w{x,t) ,U {x ,t)) . 
n— >oo 

Proof. 1. This result follows immediately from the Chapman-Kolmogorov 
equations and the fact that the fundamental solution of Ci gives the tran- 
sition densities of the Markov process X [see, e.g., Friedman (1975), The- 
orem 1.6.5.4, page 149]. From the definition of r„ and A„(6), it is easy to 
see that {xn >ti} = n and {t„ >t} = n»=i{^tfe > ^(4)} when t e {t^, t^+'^] . 
When t € [0, t^] , {r„ > t} = J] so the evaluation of pn,Wn, Un is trivial. When 
tG(tj„tj+i], nrn>t,Xt> x) = nXtl>h{ti),. . . ,Xt^>h{t'^),Xt> x), so 
using the transition probability density for the Markov process, we have 
Un{-,t) = j^(ti ) Un{y,t\)p{y,t\\ -jt) dy, from which we find the correspond- 
ing vun and pn- The first assertion thus follows. 

2. By the second part of Lemma 2.1, we have r < Tn+i < t^, and there- 
fore p < Pn+i ^ Pn, and w < Wn+i < Wn- It is clear from the definition that 
tn < tn+1 and so po = Un = Un+1 on (0,t°+i]. Now suppose Un+1 < Un on 
(0,t^+i]. Let t e {tn+i,tnXi\- Then t G {ti,ti'^^] for some j (the case t < 4 
is easier and handled similarly). Then 



Un+iix,t)= / Un+iiy,tn+i)p{y,tn+i;x,t)dy 

poo 

< / Un{y,tl+i)p{y,ti;x,t)dy 

o poo 

P{y,tn+i;x,t) / Un{z,ti)p{z,ti;y,t^^j^])dzdy 

roo POO 

< / Un{z,t^J / Piz,tn,y,tn+l)piy^t^+l^x,t)dydz 
Jb(ti,) J -oo 



INVERSE BOUNDARY CROSSING PROBLEM 



15 



= / Un{z,ti)p{z,tl;x,t)dz = Unix,t) 

Jb{ti) 

and Un+i < Un by induction on k. The proof that < C/„ < po is similar. 
The second assertion thus follows. 

3. The monotonicity of {Tn,Pn, Wn-, Un) implies the existence of the limit as 
n ^ oo. First WG show tlicit lini^^QQ r„ = r. For this, let w G be arbitrary, 
(i) If t{uj) = oo, then we have Tn{oj) = oo for all n G N so lim„_).oo Tn('^) = 
oo = t(w). (ii) Suppose t{uj) < oo. Note that (2.3) gives 

{r>t}= fl {Xs>b{s)} 
sgA(6)n(0,t) 

(3.4) 

= n n {Xs>b{s)}=f]{rn>t}. 

nGNsGA„(6)n(0,t) ngN 

Set t := lim„_i.oo ''"n.('^)- Then as r„_|_i > Tn > t for all n € N, we see that 
Tn{^) >t> t{uj) for every n € N. Consequently, u G HneNi'^" >t} = {T> t}. 
Hence, we must have t{uj) = t = limn_).oo Tn(a;). Combining the two cases, we 
obtain lim„_j,oo Tn{uj) = t{uj) for every w G J7. 

Next, we consider the limits of Wn and pn- When t = 0, we have w{-,0) = 
Wni-,0) = 1 -po{-,0) and p(0) = 1 =Pn(0). When t > 0, for each x G M, 

Wn{x,t) - W{x,t) = F{Tn > t, Xt > x) - F{t > t, Xt > x) 

= IP(r„ >t>T,Xt>x)<F{T<t<Tn). 

Thus, 

lim ||'u;n(-,i) - 'w(-,t)||Lcx,m) 

(3.5) < Ih^^ K(t)-J>(t)|= lE^ P(T„>T>t) 

n— >oo n— >oo 

= p('f|{r„>f>r}') =0. 

Finally, defining f/ := lim.„_^oo f^n we complete the proof of the proposi- 
tion. □ 

The approximating functions Un introduced in the previous proposition 
are expressed in terms of the transition densities of the diffusion X. From 
an analytic point of view, they are obtained step by step, for i = 1, 2, . . . , by 
solving the diffusion equations CiUn = in the set M x (t^,^^'''^] with initial 
values U{-,tl^+) = U{-,t) • X{b{p„),oo)i')^ where xa{x) is the indicator function 
of the set A. In the sequel, C and Ci are the differential operators introduced 
in (1.9) and (1.6), respectively. Recall the notation Qh := {{x,t)\x > b{t),t > 
0}. When b: [0, oo) — )■ [—00,00) is upper-semi-continuous, the set Qb is an 
open set with [6(0), 00) x {0} as its "initial" boundary. 
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Proposition 5. Let b: [0,oo) — )• [—00,00) be upper- semi- continuous and 
{p, w) be the survival probability and survival distribution associated with b, 
defined in (2.4)-(2.6). Then there exists a function U such that the following 
hold: 

/oo rco 
U{y,t)dy= / U{y,t)dy, 
'OO Jb(t) 

POD 

(3.7) w{x,t)= / U{y,t)dy 

J X 

(3.8) Vj;GR,t>0,0<C/</)o,0<'w<l-po, 

(3.9) ||tf;(-,t)+po(-,i)-l|lL-(R)<l-p(t) Vt>0, 

CiU Cvo<Q in Mx (0,00), 

(3.10) 

£if7 = 0, Cw = Q in Qh, 
where the inequalities in (3.10) are understood in the sense of distributions. 

Proof. Let U = lim„_s.oo Un- Since po>Un> Un+i > 0, using the Dom- 
inated Convergence theorem and the identity Wn{x,t) = Un{y,t) dy we 
obtain w{x, t) = U{y, t) dy for every x € M and t>0. Since t{uj) > t im- 
phes Xt{uj) > b{t), we see that w{x,t) = w{—oo,t) = p{t) for every x < b{t). 

It is clear that <U < po and <w <1 — po. Also, for t> 0, 

wix, t) = F{Xt > x) - P(r <t,Xt>x)> F{Xt > x) - P(t < t) 

= [l-po{x,t)]-[l-p{t)]. 

Thus, < 1 - w{x,t) -pQ{x,t) < 1 -p{t) or \\w{-,t) +Po{-,t) - 1||l=o(r) < 
l-p{t). 

It is useful to note that for each t > 0, both Wn{-,t) and [/„(•, t) are 
smooth functions. In addition, as functions of {x,t), Wn and C/„ are smooth 
in M X (0,00) \ \J^i{-oo,bnitn)] X {t^}. In particular, 

CiUn = 0, Cwn = in Qb:={{x,t)\x> bit), t>0}. 

Since both {C/nlneN and {w)n}neN are uniformly bounded in any compact 
subset of Qb, it then follows from standard results on parabolic partial differ- 
ential equations [see Friedman (1964), Theorem 3.11, page 74, and Theorem 
3.15, page 80] that uj,U £ C°°iQb) and Cw = and CiU = in Qb- 

The set of discontinuities of Un and vun is IJi^^i— 00, b{tn)] x {i^}. In 
particular, 

Wn{-,t) = Wn{-,t-), [/„(•, t) = [/„(•, t-), 

p{t)=p{t-) VtG(0,oo), 
Un{x,t+) = 0, Wn{x,t+)=pn{t) Vx < 6(i), t G {^JieN- 
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Denote by 6{- — s) the Dirac measure with mass at s and by XA the char- 
acteristic function of the set A. Then in the sense of distributions, we find 
that 

oo 

CiUn = Y,[Un{x,-)]itS{t-a 
i=0 
oo 

= -5Z^"(^'*n)'^(i-4)X(-oo,fe(tj,)](3^) <0 inMx (0,oo). 

i=0 

oo 

CWn = y^[Wn{x,-)t}tHt-tn) 

i=0 

OO 

= - Y.^Wn{x,ti) - pn{ti)]6{t - ti^}X(^_^Xt^^)]{x) <0 inMx (0,oo). 

i=0 

Sending n — ?> oo we find that Cw < and CiU < in M x (0,oo) in the 
sense of distributions. This proves (3.10) and also completes the proof of the 
proposition. □ 

4. Viscosity solutions and boundary crossing probabilities. In this sec- 
tion, we show that the survival distribution w defined in (2.6) is the unique 
viscosity solution of the time dependent Kolmogorov forward equation (1.8). 
As mentioned above, the use of viscosity solutions is necessitated by the fact 
that the boundaries arising from the solution to the variational inequality 
for the inverse boundary crossing problem do not have sufficient regularity 
for us to employ classical solutions. We do note however, that when a classi- 
cal solution exists, it gives the unique viscosity solution. Consequently, the 
classical solution of the partial differential equation (1.8), if it exists, is the 
survival distribution function associated with b, defined in (2.6). Since clas- 
sical solutions of (1.4) are obtained from classical solutions of (1.8) via the 
transformation U = —dw/dx, we also see that a classical solution of (1.4), 
if it exists, is the survival probability density of the first boundary crossing 
problem that we want to calculate. 

For simplicity, we work with b in the class Bq so that the survival proba- 
bility associated with b is continuous on [0,oo). Furthermore, we work with 
the function w which is monotone in the spatial variable, and smoother than 
U. The following definition is based on the differential inequalities/equalities 
in (3.10). 

Definition 2. Let b ^ Bq. A viscosity solution (for the survival distri- 
bution) of the boundary crossing problem associated with 6 is a function w 
defined on M x (0, oo) that has the following properties: 
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1. w£ C{R X (0,00)); limj\_,o \\w{-,t) +po{-,t) - 1\\l^(r) =0; < u; < 1; 

2. w{x,t) = w{b{t),t) Vx < b{t),t > 0; w{x,t) < w{b{t),t) Mx > b{t),t > 0; 
Cw = in Qb] 

3. If for a smooth if, point x € M and time t > 6 > 0, the function ip — w 
attains a local minimum at (x, t) on [x — (5, x + 5] x [t — d,t], then Cip{x, t) < 
0. 

We define one-sided time derivatives by 

d+(p(x,t) , (j)(x,t + At) - (p(x,t) 

K = hm , 

dt At\o At 

d-cp(x,t) , (j)(x,t) - (j)(x,t- At) 

K = hm . 

dt At\o At 

Denote by and Cf the operator C and Ci with time derivative replaced 
by the above one-sided derivative. Then from the expression of Un in (3.3), 
we see that C^Un = C~Wn = in M x Thus, we have 

(4.1) £~C/„(x,t) = 0, C'^Wn{x,t) = ^ V(x,f) gR x (0,00). 

In the uniqueness proof in the following theorem it is convenient to work 
with the special case cr = 1. This can be done without loss of generality by 
considering the transformation 

r 1 

Y(x,t):= — -dz VxGM,t>0, 

^ ' Jo <^{z,t) 

(4.2) 

Yt:=Y{Xt,t) Vt>0. 

The change of variables (x,t) — ?> (y,t) via y = Y(x,t) is smooth and invert- 
ible. Also, by Ito's lemma, 

dYt = fi{Yt,t)dt + dBt. 
Here, denoting by x = X(2/, t) the inverse of y = Y(x, t). 

Under the transformation, a boundary b for {Xt} is transformed to a bound- 
ary b:t € [0,00) Y{b{t),t) for {Yt}. Similarly, a boundary b for {Yt} is 
transferred back to b:t £ [0, 00) — )■ X(6(t), t). 

Theorem 3. Assume that b e Bq. 

1. The survival distribution associated with b defined in (2.6) is the unique 
viscosity solution for the survival distribution of the boundary crossing 
problem associated with b. Consequently, the survival probability p = V\b\ 
can be evaluated by p{-) = w{— 00, ■) where w is the viscosity solution of 
the boundary crossing problem associated with b in Definition 2. 
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2. If w is a classical (i.e., w + po e C(M x [0,oo)), dxW G C(M x (0,cxd)), 
dtw,d'^w G C{Qb)) solution of (1-8), then w is the survival distribution 
of the boundary crossing problem associated with b. If U is a classical 
(i.e., U - po€ C(M. X [0,oo)), dtU,d^U G C{Qb)) solution of (14), then 
w{x, t) := U (y, t) dt is the survival distribution of the boundary cross- 
ing problem. 

Proof. Existence. Let 6 G i?o and T,p,w be defined as in (2.4)-(2.6). 
Then p = Vlb] G C([0, oo)). We show that w is a viscosity solution in the 
sense of Definition 2. First, we show that w is continuous. Fix rr G M. For 
any t>s>0, 

w{x, t) — w{x, s) 

= P(r >t,Xt>x)- P(r >s,Xs> x) 

= P(r >t,Xt>x)- P(r >t,Xs> x) - F{t>T> s, Xs > x) 

= F{T>t,Xt> x> Xs) -¥{T>t,Xs>x> Xt) -¥{t>T>s,Xs>x). 

Note that P(t > r > s, Xs> x) <¥{t> t > s) = p{s) - p{t) so we have 

\w{x, t) -w[x,s)\< ¥{Xs >x>Xt)+ ¥{Xt >x>Xs) + \p{s) - p{t) \ . 

Since p is continuous, sending t — > s or s — )• t we conclude that w{x,-) is 
continuous in (0, oo). Next, for x <y and t > 0, 

< w{x, t) - w{y, t) = F{T>t,y>Xt> x) 

< P(y >Xt>x)= po{y, t) - po{x, t). 

Thus, w{-,t) is continuous, uniform in t G [e, oo) for any e > 0. In conclusion, 
w G C(M X (0,oo)). Recall from (3.8) that \\w{-,t) +po{-,t) - 1\\l^(m.) < 1 - 
p{t). The continuity of p on [0,oo) then implies that limt\^o||w^(') *) +Po(') *) — 
1|Il°°(R) = li™t\o(l ~ pit)) =0- Thus, w satisfies the first requirement of 
being a viscosity solution. 

Remark 4.1. The continuity of the survival probability p plays a central 
role in the proof here. In a more general case, that is, b ^ Bq, w is not 
continuous so the definition of a viscosity solution needs to be revised. To 
avoid such technicalities, we take the simple case that b & Bq. The work of 
Cheng et al. (2006) does allow discontinuous survival probabilities. 

Note that r(a;) > t implies Xt > b{t). Hence, w{x,t) = P(r >t,Xt>x) = 
w{b{t),t) when x < b{t). Also, since U = -dw/dx > 0, C/ ^ and ^if/ = 
in Qh, we have ^7 > in Qi,. In particular, if U{x,t) = 0, with {x,t) G Qb, 
then the strong maximum principle [Friedman (1964), Theorem 3.5, page 39] 
implies that U{y, t) = for all y such that {y, t) G Qb [and therefore all ?/ G M, 
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as it is easy to see that U{y,t) = for {y,t) ^ Qh]. This is a contradiction, 
since the Dominated Convergence theorem impUes that 



with the apphcation of Dominated Convergence justified by the bounds pQ > 
Un>0 from part 2 of Proposition 4. Thus, w{-,t) is strictly decreasing in 
(6(t),oo), so w{x,t) < 'w{b{t),t) for all x > b{t). Finally, from (3.10), we know 
Cw = in Qb. Thus, w satisfies the second requirement of being a viscosity 
solution. 

Remark 4.2. Recall that Vlb] = V[b*] for any boundary 6; that is, under 
our nonconventional definition of default time and survival probability, both 
the original boundary b and its upper-semi-continuous envelope b* produce 
the same crossing time, survival probability, and survival distribution. Here, 
we needed to use the upper-semi-continuous representation of the barrier so 
Qb is open and w{-,t) is strictly decreasing for x > b{t). 

We now verify the third requirement for w being a viscosity solution. 
Assume that ip — w attains a local minimum at (x, t) on A:= [x — 6,x + S\ x 
[t — 6, t] where (p is smooth and t > 5 > 0. We want to show that Cip{x, t) < 0. 
We follow a standard technique for viscosity solutions. First, we modify (p 
to a new smooth ip so — w attains a strict local minimum value zero at 
{x,t) on A. The function is defined by 



Then ip{x, t) — w{x, t) = and Cij){x, t) = Cip{x, t). That (p — w attains a local 
minimum at (x,t) implies 

'>p{y,s) -w{y,s) 



= {x- yf5-^ + {t- sf6~^ + [p{y, s) - w{y, s)] - [ip{x, t) - w{x, t)] 
>ix-y)^6-^ + {1-3^6-^ V(y,s)GA 



Thus, ip{y,s) —w{y,s) attains on yl a strict local minimum, being zero, at 



Using a standard viscosity solution technique, the differential inequality 
C(p{x,t) < is obtained by comparison of ip with smooth approximations of 
viscosity solution candidates. Here, we choose the smooth approximations 
to be {wn} introduced in Proposition 4. For each positive integer n, let Wn 
be defined as in (3.2) in Proposition 4. Then Wn is upper-semi-continuous 
on M X [0, oo), so ip — Wn attains a local minimum, on the closed set A = 
[x — 6,x + 6] X [t — 6,t]. On the parabolic boundary of A, we have ip — w >1 
and w — Wn > —I so ip — Wn > 0. At (x,t), ip — 'Wn = w{x,t) — Wn{x,t) < 0. 




i^{y, s) := ip{y, s) + {x - yfS'^ + {t - sf5~^ + [w{x, t) - ip{x, t)]. 



{x,t). 
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Hence, the minimum is attained in {x — 6,x -\- 6) x (t — 5,t]. We denote by 
{xn,tn) an arbitrary local minimizer of ip — Wn in A. 

That ip — Wn attains a local minimum at {xn,tn) implies that at {xn,tn), 
dil^/dx = dwn/dx, dtp/dt < d~Wn/dt and d'^ip/dx'^ > d'^Wn/dx'^. Hence, 

Ctp{Xn, tn) < C~Wn{Xn,tn) =0. 

In order to take the limit, we want to show that {xn,tn) — >■ {x,t) as n ^ cc. 
Intuitively this is obvious since ijj — w attains a strict local minimum at (x, t) 
and Wn^w (uniformly). Since '4>{y., s) — w{y, s) > {x — y)^6~'^ + ~ 
with {y,s) = {xn,tn), we have 

m{{Xn-x)^S-^ + {tn-tf6-^} 

n^oo 

< lim {lp{Xn,tn) -w{Xn,tn)} 
n— ^-oo 

= lim {[lp{Xn,tn) - Wn{Xn,tn)] + [Wn{Xn,tn) - w{Xn,tn)]} 



n— >oo 



< \iin < \ih(x,t) — ■Wn(x,t)] + max \wn — w\ 
" • - ' [x-S,x+S]x[t-5,t] 



= tp{x,t) — w{x,t) = 0. 
Here, we have used the uniform convergence of Wn — > w, derived as follows: 

< Wn{x, s) — w{x, s) 
= P(r„ >s,Xs>x)- P(r >s,Xs> x) 
= P(r„ >s>T,Xs>x)< P(r„, > s > r) 
= P(rn > s) - P(r > s) = pn{s) - p{s). 

Thus, we have ||i(;„(-,s) — w{-, s)\\ioo(^-j <p„(s) —p{s). Since Pn,P are con- 
tinuous, and pn \ p, the point-wise convergence of p„ — > p implies local uni- 
form convergence, that is, lim„_j.oo \\Pn —p\\l°°([o.t]) = 0- Thus, lim„_>.oo \\wn — 
w\\l^{Rx[o,t]) = for any T > 0. Hence, lim„^oo(a;„, tn) = ix,t). Finally, this 
implies Cip{x, t) = Cip{x, t) = lim„,_5.oo ^ipixn,tn) < 0. 

Uniqueness. We can assume without loss of generality that a = 1, since 
otherwise we can work with the process {lf}t>o defined in (4.2). In terms of 
our viscosity solution, it means that we make a smooth change of variable 
{x,t) {y,t) via 

r 1 

y = Y{x,t):= / -^—dz. 





In the new variables, we are working on the function w{X{y,t),t) and the 
barrier is b{X(y, t)) where x = X(y, t) is the inverse of y = Y{x, t). Retaining 
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the notation {x,t) as independent variables, we can assume that C = dt — 
\d1 + fi{x,t)dx- We denote 

^ = IIA'IIl°°(Rx[0,oo)) + l|f^a:/^||L°°(Rx[0,oo))) 
Rii) = \\Po\\L°°(Rx[t,oo)) Vt>0. 

We note that R{t) < oo by the standard Gaussian upper bound on the 
fundamental solution of Ci [see Friedman (1964), page 24]. Let w be the 
survival probability of the boundary crossing problem. Then \dxw\ < po is 
uniformly bounded in M x [tQ,oo) for any to > 0. Let w be an arbitrary 
viscosity solution of the boundary crossing problem. We want to show that 
w = w. 

Suppose w ^ w. Then there exists xq € M,to > > such that either 
w{xo,tQ) > w^XQjto) + 6(5 or ?I;(xo,io) > 'w{xo,tQ) + 65. In the former case, 
we set (101,1112) = {w, w) and in the latter case we set {wi,W2) = {w, w). Then 
both wi and W2 are viscosity solutions and 

wi{xo,to) - W2{xo,to) > 6(5 > 0. 

By spatial translation, we can assume, without loss of generality, that b{t) < 
for aU te[0,to]. 

We now fix a constant e satisfying 

0<e, eto + e^^;^ <niin((5,l), + 4e^M < e/A. 

We need another small positive constant i] determined as follows. By 
the second property of viscosity solutions, we can find ti € (0, io) such that 
\\wi{-,ti) + pQ{-,ti) — l\\ioc(^^ < 6, 1 = 1,2. Since po(' 5*1 ) is uniformly con- 
tinuous on M and W2 is continuous at {xo,to), there exist r]Q > such that 
for every 7? G (0,r?o], \w2ixo,to) - W2{xo + r],to)\ < 6 and \\po{-,ti) - poi- + 
r/, ti) ll^ooj-jg) < (5. The latter inequality implies 

\\wi{-,ti) -W2{-+7],ti)\\L^(^) 

< \\wi{-,ti) +po{-,ti) - 1||loo(k) 

+ \\w2{- + r],ti) +Po{- + r],ti) - 1||loo(ir) 

+ \\Po{-,ti) -Po(- + ??,ii)||L->(R) < 3(5. 

Now, we fix an t] € (0, ?/o] such that 

0<M??[i?(ii) + 4e2] <e/4. 

Consider the continuous function 

= wi{x,t) - W2{x + r],t) - et- e^x^, 2; G ]R,t G [ii,to]- 

Note that $(xo,to) = [wi{xQ,tQ) - W2{xQ,tQ)] + [w2{xq, to) -W2{xq + 1^,10)] - 
eto — e^a^o > 6(5 — (5 — (5 = 4(5. On the other hand, when t = ti, ^{x,ti) < 
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\\wi{-,ti) — W2{- + ?/,i2)||L°°(R) < 35. Hence, there exists (x*, i*) e M x (ii,to] 
such that $ attains at (x*,^^,) the global positive maximum of <I> on R x 

[hM- 

We consider two cases: (i) < — ??; (ii) b{t^) - r]. 

Suppose (i) X* < — r/. Then we have wi{x^,t^) = wi{b{t^),t^) and 
tf^2(a;* + ^,t*) = W2{b{t^),U). Consequently, since ||tf2(-, ||l°°(r) = 'W2{b{t^), 
t*) = W2{x* + r/, t^) and x* < < 0, we obtain 

= e'^(x^ - b{t^f) + (tt;2(x* + r/, t^) - W2ib{t^) + rj,t^))>0 

contradicting the maximality of <&(x*,t*). 

(ii) X* > — r/. Set ip{x, t) = 'W2{x + r], t) + et + e^x^. Then ip — wi = — $ 
attains at (x*,t*) a minimum over M x [toi^i]- Since x^, > 6(t*) — r/, we see 
that W2{- +??,•) is smooth in a neighborhood of (x*,t^,). Then if is smooth 
in a small neighborhood of (x*,t^,) and (p — wi attains a local minimum at 
(x*,t*). Since wi is a viscosity solution, we must have Cip{xt,,t^) < 0. 

Now, we calculate Cip{x^,t:t). Using Cw2 = in and the fact that a is 
assumed to be a constant, we have 

>C(/?(x*,i*) = e — + 2e^x*/i(x*,t^,) 

+ [/i(x*,t^,) -^(x* + ??, i*)]92,.W2(x* + r/,t^). 

First, we note that $(x^,, t*) > <I>(xo,to)i so et* + e'^x^ < 2 + eto + e^Xo < 3. 
This implies that e^|x*| < 2. Hence, 

Cip{x^,t^) >e - - Ae^M - Mr]\dxW2{x* + ??, 

To estimate dxW2{x.f + r],t^), we consider two situations. 

(a) W2 = w is the survival distribution function. Then [[^^^'"^(a^*, i*) || < 
Po(a^*, t*) < R{ti)- 

(b) W2 = w. Then ^2 is differentiable at (x^=,t*) and is Lipschitz 
continuous with Lipschitz constant ||C^(-,i*)||L°°{R)- Hence, sending /i \ in 
[($(x* ± - ^>(x*,t*)]//i > we derive 

\dxW2ix^,t^)\ < ||c?a:t(;(-,t*)||iao(iR) + 2e^|x*| < + 4e^. 

Thus, in either case, we have 

Cip{x^,Q > e - - Ae'^M - Mr]{R{ti) + Ae'^} > e/4 > 

by our careful choices of e and r/. This contradicts C(p{x^,,t^) < 0. The con- 
tradiction implies that we must have wi = W2- Thus, the viscosity solution 
of the boundary crossing problem is unique. 

Proof of the second assertion. The equivalence of classical solutions of 
(1.8) and (1.4) via U = —dw/dx is trivial. Here we show that a classical 
solution of (1.8) is a viscosity solution. 
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Assume that u; is a classical solution of (1.8). Then U = —dw/dx is a 
classical solution of (1.4). Applying the maximum principle to U and po ~ 
U on Qb, we find that < C/ < po- Also, the strong maximum principle 
[Friedman (1964), Theorem 3.5, page 39] shows ?7 > in [if U (x, t) = for 
(x, t) G Qfc, then U {y, s) = for all (y, s) in Qi, with s <t, contradicting the 
initial condition at time 0]. Hence, w is monotone in x and w{x, t) < w{b{t),t) 
for all x > b{t),t > 0. In addition, for each tQ > 0, comparing U with the 
solution of CV = in M x (tg, oo) with initial value V{-,to) = U (•, to) we see 
that U<V onRx [to, oo) so U (y, t) dy < V{y, t) dy = /jj U{y, to) dy for 
every t > to- This implies that p{t) := 'w{—oo,t) = w{b{t),t) is a decreasing 
function of t. 

Next, since w + po £ C(M x [0,oo)), dxW < 0, dxPo < 0, and u'(oo,0) = 
and w^—oo, 0) = 1, one can show that lim^x^o ll^(') *) +Po('i — 1||l°°{r) = 0. 
Thus, w satisfies the first and second requirements of a viscosity solution in 
Definition 2. 

To verify the third requirement in Definition 2, suppose (p is smooth, 
x (^M,t > 6 > and ip — w attains a local minimum at (x, t) on [x — 6,x + 
6] X [t — 6,t]. We want to show that C<p{x,t) < 0. We consider two cases: (i) 
X > b{t) and (u) x < b{t). 

(i) Suppose X > b{t). Then (x,t) is an interior point of Qb, in which w 
is smooth. Since ip — w attains a local minimum at (x, t) on [x — 5, x + (^] x 
[t — 5,t], we have dtip{x,t) < dtw{x,t),dx(p{x,t) = dxw{x,t) and dx(p{x,t) > 
dxw{x,t). This implies that Cp{x,t) <Cw{x,t) = 0. 

(ii) Suppose x < b{t). Then dxip{x,t) = dxw{x,t) = 0. Note that 99(x — 
z, t) — (p{x, t) > w{x — z,t) — w{x, t) =0 for all z>0. This implies, since ip is 
a smooth and dxip{x,t) = 0, that dxxP{x,t) > 0. 

To complete the proof that Cip{x,t) < 0, it suffices to show that dt^p < 0. 
Suppose on the contrary that dtip{x,t) > 0. Then there exists e € {0,S) such 
that ip{x, t — s)< ip{x, t) for all s £ {0,£]. As p — w attains a local minimum at 
{x,t), we see that ■w{x,t — s) < w{x, t) — ip{x, t) + ip{x,t — s) < w{x, t) =p{t) < 
p{t — s) for all s G (0, e] . Thus, [x, oo) x [t — e, t) C Qh- Since w is monotone, 
we also have w{y,t — s) < p{t) = w{x,t) for all y > x,s € [0,e]. That is, w 
attains at (x,t) a local maximum over the region [x,oo) x [t — e,t]. Hence, 
applying Hopf's lemma [Protter and Weinberger (1967), Theorem 3.3] for w 
on [x, oo) X {t — s,t), we have Wx{x, t—) < 0, which contradicts the definition 
of a classical solution that dxw{x,t—) = dxw{x,t) = 0. Thus, we must have 
dtp{x,t) < 0. Together with dxp{x,t) = 0,dxxP{x,t) > 0, we conclude that 
Cp{x,t)<0. 

Hence, w is a viscosity solution. Applying the conclusion of the first asser- 
tion, we then see that w is the survival distribution of the boundary crossing 
problem associated with b. □ 
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Remark 4.3. (i) The addition of the term e^x^ confines our attention 
in searching for the maximum of $ to a compact set [— 2e~'^, 2e~^] x [0, io]- 
So we indeed only need w +po G C(M x [0, oo)) and w + po — 1 = on M x {0} 
in the definition of viscosity sohitions. 

(ii) Fix ^2 > ^1 ^ 0. In the proof, if we take W2 to be a solution of Cw2 = 
in M X {ti,t2) subject to initial condition W2{-,ti) > Then following 

the proof we see that supigx[ti,j2)(''^ ~ ''^2) > is impossible. Thus, we have 
w <W2 on R X [ti, ^2)- This is a simple version of the comparison principle in 
the theory of viscosity solutions. This result will be used in the next section. 

5. Viscosity solutions for the inverse boundary crossing problem. The 

inverse boundary crossing problem is to find 6, for a given p, such that p is the 
survival probability associated with b. In this section, we prove that for any 
p (z Pq, from the viscosity solution of the variational inequality (1.11), we can 
find an unique b € Bq such that the resulting p gives the survival distribution 
of the first time that X crosses b. Since the forward problem maps Bq to Pq, 
we study, for simplicity, the inverse problem for p Pq, though in Cheng et 
al. (2006) the variational inequality (1.11) was considered for more general 
survival functions. 

5.1. Viscosity solutions. In general, classical solutions of the variational 
inequality (1.11) for the inverse problem may not exist. In Cheng et al. 
(2006), viscosity solutions were introduced, and it was shown that for any p 
satisfying 

(5.1) p{0+) = l>p{s)=p{s-)>p{t)>0 Vt>s>0, 

there exists a unique viscosity solution. Prom this solution, we can define a 
boundary b such that Qb = {w <p}, and consider it as a candidate for the 
solution to the inverse boundary crossing problem. To verify that b is indeed 
a solution, we show that t;; is a viscosity solution to the direct problem 
(1.8), and then appeal to Theorem 3 to see that w and p give the survival 
distribution of the first time that X crosses b. 

When p Pq, we know a priori that the unique solution of the variational 
inequality is continuous so many technicalities in Cheng et al. (2006) re- 
garding the definition, existence, and uniqueness of viscosity solutions can 
be ignored. In particular, the viscosity solution introduced in Cheng et al. 
(2006) can be reformulated (removing those specifics that take care of dis- 
continuities) as follows. 

Definition 3. Let p G Pq be given. A viscosity solution for the survival 
distribution of the inverse boundary crossing problem associated with p is a 
function w defined on M x (0, 00) that has the following properties: 
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1. w+poe C{R X (0,00)), limt^o \\Po{-,t) + w{-,t) - 1||lcx,(ir) =0; 

2. < w{x, t) < p{t), V(x, i) G M X (0, 00) and £w{x, t) = in the set Q := 
{ix,t)\t>0,wix,t)<pit)}; 

3. if X € M and t> 6 > 0, and is a smooth function such that ip — w attains 
at {x, t) a local minimum on [x — 5, 2; + (^] x \t — 6, t], then Cip{x, t) < 0. 

A viscosity solution of the inverse boundary crossing problem associated 
with p is the function b given by 

(5.2) b{t):=mi{xeR\w{x,t)<p{t)} Vt > 0, 

where it; is a viscosity solution for the survival distribution of the inverse 
boundary crossing problem associated with p. If there is a unique viscosity 
solution, we denote b = B[p]. 

The remainder of this section is devoted to a proof of the main result of 
the paper, Theorem 1, stated in the Introduction. 

Proof of Theorem 1. The fourth assertion follows from the first as- 
sertion and the following facts which are easy to verify: a classical solution 
of (1.11) is automatically a viscosity solution, and if {U,b) is a classical so- 
lution of (1.4), then {w,b) with w defined by w{x,t) = U{y,t) dy is a 
classical solution of (1.11). We divide the proof of the first three assertions 
into several parts. 

Existence and uniqueness of a viscosity solution. The proof of the exis- 
tence of a unique viscosity solution, together with the properties (1.14), is 
the major result of Cheng et al. (2006) and hence is omitted here.^ It is 
important to note that, by the monotonicity of w in the spatial variable and 
the definition of b in (5.2), we have 

Q := {{x,t)\t > 0,w{x,t) <p{t)} = Qb := {{x,t)\t >0,x> b{t)}. 

Weak regularity of the free boundary. The regularity of the free boundary 
b = B[p\ defined by (5.2) was not discussed in Cheng et al. (2006). Here, 
under the assumption that p € Pqi we establish a very basic regularity result 
on b. In particular, we show that b Bq. 

We begin by showing that b{t) < 00 for every t > 0. Indeed, if b{t) = 
00, then by the definition b{t) = mi{x\w{x,t) < p{t)} we see that w{x,t) = 
p{t) for all X € M. Since \\mx^oow{x,t) = (recalling u; < 1 — po)j we see 



^In Cheng et al. (2006), it was assumed that P(Xo = 0) = 1. The same techniques can be 
applied to prove existence and uniqueness for more general initial distributions considered 
here. 
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that w{-,t) = 0. This contradicts the assumption that p € Pq, since p € Pq 
guarantees p{s) > for all s € [0, oo). 

Next, we show that Xq > 6*(0) = limj\^o b{t) almost surely. To see this, we 
use the estimate w <1— po to derive p{t) = w{b(t),t) < 1 — pQ{b{t),t). This 
implies that limt\oPo {b{t),t) < limf_j.o(l — = since p € Pq gives p € 
C([0, oo)) and p{0) = 1. Now suppose to the contrary that P(Xo < 6*(0)) > 0. 
Then there exists > and e > such that po {h* (0) - 25, 0) =¥{Xq< b* (0) - 
25) > 2e. Consequently, there exists to>0 such that po{b*{0) -6,t) = F{Xt < 
b* (0) — 6)>e for aU t £ [0, to] • However, by the definition of b* (0) , there exists 
a sequence tk\0 such that lim/c-yoo b{tk) = b*{0). For all sufficiently large k, 
wehave6(tfc) > 6* (0) -(^ which implies that po(&(*A:), ifc) >p{b*{0)-d,tk) > e, 
so we obtain limt\oPo{b{t),t) >e. This contradicts limt\opo{b{t),t) = 0. 
Hence, we must have P(Xo < 6*(0)) = 0, that is, Xq > b*{0) a.s. 

The next step is to show that b := B[p] is upper-semi-continuous on [0, oo). 
First of all, the definition 6(0) := limsup^\^o ^(^) implies that b is upper-semi- 
continuous at t = 0. Next, let t > be arbitrary. We consider two cases: (i) 
b{t) > -oo, (ii) b{t) = -oo. 

(i) Suppose b{t) > -oo. Fix any e > 0. Then p{t) - w{b{t) +e,t)> 0. By 
continuity, p{s) — w{b{t) -|- e, s) > for all s in a neighborhood of t. This 
implies that b{s) < b{t) + e for all s in a neighborhood of t. Consequently, 
\\m.s^t b{s) < b{t) -|- e. As e > is arbitrary, we have lims_j.f b{s) < b{t). 

(ii) Suppose b{t) = —oo. Then for any M > 0, we have p{t) — w{—M, t) > 0. 
Consequently, p{s) — w{—M, s) > for all s sufficiently close to t. Hence, 
b{s) < —M for all s sufficiently close to t. This implies that lims^tb{s) < 
—M. As M can be made arbitrarily large, we hence see that Yimg^t b{s) = 
-oo = b{t). 

In conclusion, b: [0,oo) [—00,00) is upper-semi-continuous. 

Let Q := {w <p} := {{x,t) e Mx (O,oo)|'u;(x,t) <p(t)} and Qb := {{x,t) G 
Mx (0, oo)|j; > b{t)}. Then Q = Qb and U := -dxW > in Qb- Indeed, Q = Qb 
follows from the definition oi b = B\p] in (5.2) and the monotonicity of w in 
the spatial variable [see Cheng et al. (2006)]. In addition, since b is upper- 
semi-continuous and bounded above in any compact interval in [0, 00), any 
points {xi,ti), {x2,t2), ti < t2 in Qb can be connected by a smooth curve 
X = h{t) in Qb defined on t € [ti, i2] such that h{ti) = Xj. Applying the strong 
maximum principle [Friedman (1964), Theorem 3.5, page 39] to U := —dxW 
{U >0 hy the monotonicity of w), we conclude that U > in Q = Qb [an 
elementary argument shows there cannot exist at2> such that [/(•, ti) = 
in Qb for all ti < 12]. 

Next, we show that b = b* = b*_. Upper-semi-continuity (6 = 6*) was shown 
above, so it remains to prove that 6(t) = lim^/if 6(s) =: 61 (t) for every t > 
0. Let t > 0, and suppose b{t) ^ b*_{t). Then b{t) = b*{t) > b*_{t). Set 6 = 
[b{t) - b*_{t)]/4:. By the definition of b*_{t), we can find e > such that 6(s) < 
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61 (t) + 6 for all s£[t-2e,t). Then D:=[b*^ (t) + 6, b{t)] x[t-2e,t)isa subset 
of Q := {w < p}. Since we know that CiU = and U = —dw/dx > in Q. 
We can apply the Harnack inequality on the cube (61(t) + 5, x {t — 2e,t) 
to conclude that there exists a positive constant tj > such that C/ > 77 in 
[b*_ (t) + 2(5, b{t) -6]x {t-e,t). Consequently, 



Sending s /^t we then conclude that w{b* + 26, t) > w{b{t) — 6,t) + 67] > 
p{t) + which violates the requirement that w <p for a viscosity solution. 
Hence, we must have b{t) = b*_{t). In summary, b = b* = b*_. This also implies 
that 



Finally, to show that b & Bq, it remains to show that the survival prob- 
ability p := Vm associated with b has the property limj\^oP(0 — 1- For 
this, we consider the sequence {m„}, associated with 6, defined in Proposi- 
tion 4. It follows from a (viscosity solution) comparison principle, applied 
iteratively to M x (t^,*^"*"^] {t^ '■= 0) for i = 0,1, . . . , that Wn > w; see Re- 
mark 4.3. Taking the limit, we find that w < lim„_;,oo "w^n • This implies that 
p{t) = w{—oo,t) < lim„_s.oo OO) ^) = Pit)- Since our assumption p € Pq 
implies that limf\^Qp(t) = 1, we also know that lim(\^oP(0 — 1- Thus, we 
have shown that 6 € i?o • 

Verification that the boundary derived from the variational inequality has 
the required crossing time distribution. Given p £ Pq, let b = B[p] be the 
boundary derived from the unique solution of the variational inequality 
(1.11). We need to show that V[b] = p, that is, that b is truly a solution 
of the inverse boundary crossing problem. Summarizing, this means that we 
want to show that {V o B) [p] = p for every p G Pq. 

Let w be the unique viscosity solution for survival distribution of the 
inverse problem associated with p as given in Definition 3. Define b = B[p] as 
in (5.2). Let p = V[b]. We want to show that p = p.lt is enough to show that 
w is a, viscosity solution of the survival probability for the boundary crossing 
problem associated with b, since in this case part 1. of Theorem 3 yields that 
p{t) = w{—oo,t), while taking limits as x goes to —00 in (1.14) gives that 
p{t) = w{—oo,t). By checking the Definitions 3 and 2 of viscosity solutions, 
one readily sees that w being a viscosity solution in the sense of Definition 
3 implies that w is indeed the viscosity solution in the sense of Definition 



w{b*_it) + 26,s)-w{b{t)-6,s 




6rj 



VsG {t-e,t). 



= lim6(s) = lim6(s) > lim6(s). 

s~>t s /-t s\t 
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2, provided that Q := {w < p} = Qb ■= {{x,t)\x > b{t),t > 0}. But this last 
property is immediate from the monotonicity of w. We thus conclude that w 
is indeed the viscosity solution of the survival distribution of the boundary 
crossing problem associated with 5. Consequently, V[b]{t) = w^—oOjt) =p(t), 
so we have V[b] = p and p = V[b] = {V o B) [p] . 

Uniqueness of the solution of the inverse boundary crossing problem in 
the class Bq. For a given p € Pq, we have shown that is a solution of 
the original inverse boundary crossing problem. Here we show that there 
is indeed only one such b in the class Bq. To show this, it suffices to show 
that {B o 'P)[b] = b for every b € Bq, since this implies that if Vlb] = p then 
b={Bo^V)[b] = B[p]. 

Let b Bq. Define {t,p,w) as in (2.4)-(2.6). That is, p = Vlb] and w 
are the survival probability and distribution of the boundary crossing prob- 
lem associated with b. Since b is upper-semi-continuous we can derive from 
Proposition 5 and the strong maximum principle that U := —dxW > in 
(see the proof of Theorem 3). This implies that C {w <p}- Also, since 
w{x,t) = P(t '>t,Xt> x) we know that w{x,t) = p{t) when x < b{t). Thus, 

Q'b = {w<p}. 

By Theorem 3, t/; is a viscosity solution of the survival probability distribu- 
tion of the boundary crossing problem in the sense of Definition 2, associated 
with b. By checking the definition of a viscosity solution of the variational 
inequality associated with p (Definition 3), we find that w is indeed a vis- 
cosity solution associated with p = V[b]. Now, according to the definition 
of B[p] in (5.2), B[p]{t) = mi{x\w{x,t) <p{t)}. Since Qi = {w < p}, we see 
that B\p] = b. Thus b = B[p] = {B o V)[b] for every b G Bq. This completes 
the proof of Theorem 1. □ 

5.2. Continuity of the free boundary in the inverse boundary crossing prob- 
lem. In this subsection, we investigate the continuity of the free boundary 
5 = B[p] for the inverse boundary crossing problem for p € Pq. We already 
know that b is upper-semi-continuous, and since b = b*_, it cannot "jump 
up." For b to be continuous, we need to prevent it from "jumping down." 
Note the fact that if p is a constant in an open interval, then b = — oo in 
that interval. Hence, to eliminate steep drops of b we require a lower bound 
on the rate of decrease of p. We consider the following: 

The following proposition gives a sufficient condition for the boundary to be 
continuous, in the case that X is a standard Brownian motion. 
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Proposition 6. Suppose that X is a standard Brownian motion, that 
is, /U = 0, a = 1, and po{x,0) = X[o,cx,)ix)- Let p G Pq and b = B[p]. 

1. // L(p, Ti, T2) > for some positive Ti, T2 with Ti < T2, then b = B[p\ 
is continuous on (Ti,T2). 

2. Assume that L{p,0,T) > for every T>0. Then b e C{[0,oo)). 



Proof. 1. Let ti > be arbitrary. Define 

, W(b(tl) +X,t + tl) 7,^ ■,, N ,/ N 

w{x,t)= ^ ^ \' -, b{t) = biti + t)-biti), 

Pit) 



Then {w,b) is the solution of the inverse problem with initial value w{-,0) 
and survival probability p. This statement follows by an immediate appli- 
cation of the definitions. Note that w{x,t) =¥{Xt-^-ti > x,t >t + ti\T >ti). 
p{t) = P(r > t + tijr > ti). The conditional probabilities and the boundary 
from time ti on are the same as the solution of the inverse problem started 
with the initial position equal to the conditional distribution of X given that 
r > ti. 

Now let be the solution of the inverse problem with initial data 

X(-oo,o) ^iid survival probability p. Note that t/^(x,0) = 1 = w{x,Q) for x < 
and w_{x,{)) = < w{x^Qi) for x > 0. Hence, w;(-,0) < w{-,d). It then follows 
from a comparison principle [cf. the proof of Lemma 4.2 in Cheng et al. 
(2006)] that w<w and that b<b. Again, this is obvious from the probabilis- 
tic interpretation of the problem. The boundary b is the one that produces 
the hitting distribution p when the process starts at b{ti) at time ti. The 
boundary b produces the same hitting distribution with the process started 
at the conditional distribution of Xt^ given that t >ti. Since in this case 
we must have Xt^ > b{ti), we have that the boundaries b and b produce the 
same hitting distribution for the process X, with b arising from X starting 
at a higher point with probability 1. Therefore, we must have b>b. 

Thus, for 0<t< 1/2, 



6(ti + t)- biti) = b{t)>b{t) >-[! + 0{l)y -2t\og[l - p{t)] 

by the estimate on line 16, page 867 of Cheng et al. (2006). Upon noting 
that 

|log[l -p{t)]\ = |logb(ii +t) -p{ti)] - logp(ii)| 

= \\og{tp{ti + e))\ + o{i) 

= |logt|+0(l) = 0(l)|logt|. 
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where 

p{t) := limsup ~ ^^^^ € [0, oo] 

s/-t t — S 

we find that there exists a constant C{ti) such that 

b{t + h) - b{ti) > -C{ti)y^t\logt\ yt G (0, 1/2]. 

Now pick any t G (Ti,T2). Let {ti}^^ be a sequence in [Ti,t) such that 
linii^oo ti = t and linii^oo b{ti) = b{t) [recalhng b{t) = b*i{t) := limg^j 6(s)]. 
Then setting hi = [t — ti]/2, we have 

b{ti + h) > b{ti) - cyJ\hi\oghi\ V/i G [/li , 3/ii] . 

This implies that 

inf b{s) > b{ti) - cyJhi\\og hi\, 

se[t-hi,t+hi] 

SO that 

lim b(s) > b{ti) — c^hi\\og hi\. 

s—>-t 

Sending i — )■ oo we then obtain lim ,_^^ b(s) > b{t). Thus, b is lower-semi- 
continuous in (Ti,r2). Since b is also upper-semi-continuous, we see that b 
is continuous in {Ti,T2). 

2. By the first assertion, we know that b is continuous in (0, oo). At t = 0, 
since po{x,0) > for all x > 0, the proof in Cheng et al. [(2006), Lemma 4.5, 
page 865] implies that there exists a positive constant C that depends on 
L(p, 0,1/2) such that 

b{t) > -C^/t\logt\ Vt G [0, 1 /2] . 

This implies that lim ^\^Q b(t) > 0. 

We recall that 1 — w{x,t) — pQ{x,t) > 0. Evaluating this inequality at 
X = b{t) gives po{b{t),t) < 1 — p{t) for all i > 0. Since pq{x^Qi) > 0, sending 
t\0 we conclude that lim^x^o b{t) < 0. Thus, 6(0) := lim^x^o ^(0 = 0- This 
completes the proof. □ 

Note that if the hitting time density —p is everywhere strictly positive, 
then we obtain that b = B[p] is continuous. In particular, this criterion is 
satisfied by the boundary arising from the exponential distribution and hence 
provides the solution to the inverse boundary crossing problem as originally 
proposed by A. N. Shiryayev. 
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